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We have calculated the kinematical correlations between charged leptons from semilep- 
tonic decays of open charm/bottom as well as leptons produced in the Drell-Yan mech- 
anism in proton-proton collisions at BNL RHIC. Presented production rates for open 
charm and bottom were estimated by a detailed analysis of so-called nonphotonic elec- 
trons. The distributions for charm and bottom quarks pairs are calculated in the kt- 
factorization approach with the help of the Kwiecihski unintegrated parton distributions. 
The hadronization of heavy quarks is done by means of Peterson et al. fragmentation 
function. The semileptonic decay functions are found by fitting recent semileptonic data 
obtained by the CLEO and BABAR. We get good description of recent PHENIX data. 
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1. Introduction 

Recently the PHENIX collaboration has measured dilepton invariant mass spectrum 
from to 8 GeV in proton-proton collisions at y/s = 200 Gew4 It is believed that the 
main contribution to the dielectron continuum comes from so-called nonphotonic 
electrons which are produced in semileptonic decays of open charm and bottom 
mesons. To date, productions of charm and bottom was usually studied by standard 
measurements^ and pQCD calculations^ of single leptons inclusive distributions. 
Such predictions give rather good description of the experimental data, however 
the theoretical uncertainties are quite large what makes the situation somewhat 
clouded and prevents definite conclusions. 

Better statistics at present colliders gives a new possibility to study not only 
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inclusive distributions but also correlations between outgoing particles. Kinemati- 
cal correlations constitute an alternative method to study open charm and bottom 
production. It gives also a great possibility to separate charm and bottom contribu- 
tions what has a crucial meaning for understanding the character of heavy quarks 
interactions with the matter created in high energy nuclear collisions^. 

The process of nonphotonic electron production has three subsequent stages. 
The whole procedure can be written in the following schematic way: 

da e do® 
dyd z p dyd z p 

where the symbol ® denotes a generic convolution. The first term is responsible 
for production of heavy quarks/antiquarks. Next step is the process of formation of 
heavy mesons and the last ingredient describe semileptonic decays of heavy mesons 
to electrons/positrons. 



2. Formalism 

The inclusive production of heavy quark/antiquark pairs can be calculated in the 
framework of the fc t -factorization-2J. In this approach transverse momenta of initial 
partons are included and the emission of gluons is encoded in so-called unintc- 
grated gluon (parton) distributions. In the leading-order approximation within the 
fc t -factorization approach the differential cross section for the QQ or Drell-Yan pro- 
cess can be written as: 

do_ f d 2 Ki it d 2 K 2 ,t 1 |^ p ^ 

dy 1 dp lt dy 2 dp 2 td(t) 4~. ' J ir ir 1Qtt 2 (x 1 x 2 s) 2 y 

S 2 (k m + K 2 ,t ~Pi,t ~P2,t) Fifai, «i,t) Fj{%2, «!,t) , 

where J~i(xi, k\ t ) and J-j{x 2 , k 2 t ) are the unintegrated gluon (parton) distribution 
functions (UPDFs). 

There are two types of the LO 2 — !► 2 subprocesses which contribute to heavy 
quarks production, gg — > QQ and qq — > QQ. The first mechanism dominates at 
large energies and the second one near the threshold. At relatively low RHIC ener- 
gies rather intermediate x-values become relevant so the Kwiecinski UPDFs seem 
applicable in this case^. 

The hadronization of heavy quarks is usually done with the help of fragmen- 
tation functions. The inclusive distributions of hadrons can be obtained through a 
convolution of inclusive distributions of heavy quarks/antiquarks and Q — > h frag- 
mentation functions. The Peterson fragmentation functions are often used in this 
contextP. 

Recently the CLEO and BABAR collaborations have measured very precisely 
the spectrum of electrons/positrons coming from the weak decays of D and B 
mesons, respectively". These functions can in principle be calculated. This intro- 
duces, however, some model uncertainties and requires inclusion of all final state 
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channels explicitly. An alternative is to use proper experimental input which af- 
ter renormalizing to experimental branching fractions can be use to generate elec- 
trons/positrons in a Monte Carlo approach. 

The kt -factorization method is very useful to study correlations between pro- 
duced leptons because it can reproduce distributions in kinematical variables re- 
lated to the transverse momentum of initial partons. In this part of our calculations 
we take under consideration not only leptons from open charm/bottom decays but 
also leptons produced in Drell-Yan proccess, as well as leptons coming from elastic 
and inelastic processes initiated by photon-photon fusion. In the case of elastic re- 
action we follow exact momentum space calculations with 4-body phase space and 
for inelastic scattering we have applied unique photon distributions in the nucleon 
MRST200J3. 

3. Results and Conclusions 
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Fig. 1. Dilcpton invariant mass spectrum with PHENIX data and our predictions (left panel) 
and azimuthal angle correlations between nonphotonic and Drell-Yan leptons in different invariant 
mass regions (right panel). 

We have calculated dilepton mass distribution as well as azimuthal angle cor- 
relations between outgoing leptons. We get a rather good description of the recent 
PHENIX data. 
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